High-T c Superconductors in Applied Magnetic Fields Parallel to the CuO Planes: 
First Order Transition with Slow Onset of Resistivity 
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The three dimensional uniformly frustrated XY model is used as a model of a high temperature 
superconductor in an applied magnetic field parallel to the CuO-planes. Through Monte Carlo 
simulations with anisotropy n 2 = 10 on large lattices we find evidence for a first order transition. 
Earlier simulations and theoretical treatments are discussed and the experimentally found smooth 
onset of resistivity is suggested to be due to a large potential barrier against vortex line motion 
above T c present for perfect alignment of the applied field. 
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The behavior of cuprate superconductors in applied 
magnetic fields has been, and continue to be, a very ac- 
tive area of research. The most studied geometry is with 
the applied field perpendicular to the CuO-planes. In 
that case, for clean samples, a first order transition as- 
sociated with the melting of the Abrikosov lattice is by 
now well established both through calorimetric measure- 
ments jl]] and measurements of the jump in the magnetic 
induction g|. Also, a number of transport measurements 
are consistent with a first order transition. 

The geometry with the applied field parallel to the 
planes, has not been that thoroughly examined. It has, 
however, been found that the sharp increase in resis- 
tivity associated with the first order melting transition 
vanishes for perfect alignment of the field, and that the 
temperature-dependence of the resistivity is consistent 
with a continuous transition j^] . This remarkable finding 
naturally rises the question of the possibility of a contin- 
uous vortex line lattice melting, and the commonly ac- 
cepted theoretical picture is that the system undergoes a 
continuous transition analogous to the nematic to smec- 
tic transformation in liquid crystals [Q. A more recent 
experiment |Bj suggests that the transition from a vortex 
liquid to a vortex solid phase takes place in two steps, 
with an intermediate 'possible smectic' phase. 

In the first Monte Carlo study of this model || the 
authors examined a model with anisotropy rj 2 = 100 (see 
below), and argue that the behavior of both the helic- 
ity modulus and the specific heat suggest a continuous 
transition. However, a comparison between their Fig. 2 
and corresponding data for a 2D XY model of the same 
size shows a remarkable similarity, suggesting that it is 
the individual layers that dominate the behavior of this 
three-dimensional system. This makes it difficult to draw 
any conclusions based on these results. 

In this Letter we report on a study of a similar system, 
but with anisotropy rj 2 = 10. The smaller value for the 
anisotropy is chosen with the hope to get a more clearly 
visible signal from the transition. The behavior of the 
system is clearly dependent on the system size. Only for 
very large systems does it become clear that the melt- 



ing transition is first order instead of continuous. We 
also discuss problems with the suggested smectic phase 
and suggest a simple mechanism behind the apparently 
continuous onset of resistivity at the transition. The sug- 
gested scenario is directly open for experimental verifica- 
tion. 

For the study of superconductors in applied magnetic 
fields we make use of the uniformly frustrated anisotropic 
3D XY model, with the Hamiltonian 



H 



A, 
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where the sums are over the lattice points i and jl = x, 
y, and z. To include a magnetic field in the y-direction 



we choose A- 
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where D is the discrete difference operator. We use equal 
couplings in the x-y planes, J x = J y 



J xy , but a differ- 




ent coupling, J z , between the planes. The anisotropy is 
given by the parameter 



(2) 



where the relation to the physical parameters, the pen- 
etration lengths A M in the different directions, the inter- 
plane distance d and the bare vortex core size in the 
planes is given by the second equality. A derivation 
and justification of Eqs. (0) and (|^) may be found in Ref. 

To examine the behavior of this system we have sim- 
ulated Eq. (|l|) with periodic boundary conditions, / = 
1/24, and several values for the anisotropy and system 
sizes. The results reported here are for rj 2 = 10 and 
L x = L z = 48, and L\\ = 48, 96, 192, and 384. We mea- 
sure several different quantities. The superconducting 
coherence in the different directions as examined through 



the helicity moduli T X: Y z , and Tn 



To study 



the melting of the vortex lattice we measure the structure 
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factor S^kjJ where is a wave vector perpendicular to 
the field and monitor AS* = |5(K) - S(R X [K])\ § where 
K is a reciprocal lattice vector of the ordered vortex lat- 
tice, and R x reflects K through the x axis, and we average 
over the three smallest non-zero values of K. Our Monte 
Carlo simulations consist, for the larger systems, typi- 
cally of 10 5 passes through the lattice for equilibration 
followed by 1 x 10 6 through 4 x 10 6 passes for calculating 
averages. The runs on the smaller systems were typically 
somewhat shorter. 

For the geometry with the field perpendicular to the 
CuO layers it has been found that the vortex lattice melt- 
ing transition is best examined by starting with a disor- 
dered system and slowly lowering the temperature until 
the system finds an ordered configuration with a vortex 
lattice ||. However it turns out that the choice of Ln 
- the extension of the system in the field direction - is 
crucial in order to succeed with this cooling into an or- 
dered state and obtain the correct transition behavior. 
With a too large value of Ln it becomes very difficult - 
or due to limited computer resources even impossible - 
to cool into a lattice. With a too small Ln the simula- 
tions incorrectly give at hand that there are two separate 
transitions with the vortex lattice melting taking place 
at the lower temperature followed by the vanishing of T|| 
at a higher temperature. The proper choice of the length 
L|| depends on the parameters of the model, / and r\. In 
the first simulations with clear evidence of a first order 
transition, with rj 2 = 10 and / = 1/25, this length was 
chosen as L\\ — 40 M, but for isotropic systems, rf = 1 
and / = 1/20, the evidence of a sharp transition was only 
found for L|| = 128 g or L|| = 120 0. 

Guided by the experience with fields perpendicular to 
the CuO planes we therefore first tried slowly cooling 
the system, but we never succeeded cooling into a lat- 
tice. We therefore instead started the simulations with 
an artificially prepared initial configuration with a per- 
fect vortex lattice and slowly increased the temperature. 
In this way we obtained results similar to Fig. 1 of Rcf. 
M. The vortex lattice melting as probed by A5*(K) is 
shown together with T|| in Fig. |l| for L\\ = 48, 96, 192, 
and 384. 
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FIG. 1. The melting of the ground state vortex lattice for 
L x = L z — 48 and several different values for Lu. Shown is 
the structure factor together with Ty. Only for Lu > 192 
does it become clear that there is a single transition with 
simultaneous vanishing of both AS(K) and Ty. 

This data shows that one needs at least L|| ss 192 for 
obtaining evidence that the vanishing of goes together 
with the melting of the vortex lattice. The similarity with 
data from Ref. ^ for the field perpendicular to the planes 
strongly suggests a first order transition also for this case 
with the field parallel to the planes. 

From these simulations for rj 2 = 10 we can of course not 
immediately draw any conclusions regarding the charac- 
ter of the transition for larger r\. Still, we will argue below 
that both the simulation results presented in Ref. || and 
the experimentally found smooth onset of the resistivity 
H are consistent with first order transitions. Since we 
also find problems with the continuous melting scenario, 
we suggest that the transition actually is first order for 
all finite values of the anisotropy. 

Turning to the simulation results for r/ 2 — 100 ||, we 
first note that the choice L|| = 40, or 80, in the light of 
Fig. [IJ not is large enough for giving the correct behavior. 
But it also seems that their results are dominated by the 
behavior of the (weakly coupled) two-dimensional layers 
in the x-y plane. To understand the relation to the 2D 
XY model we note that it is possible to choose a gauge 
with Ai x — Ai y — 0, and the frustration included in the 
Ai z . Considered in this way, the system consists of some 
stacked unfrustrated 2D XY planes with a weak inter- 
plane coupling that at some points between the planes is 
ferromagnetic (Ai Z = 0), at some points is antiferromag- 
netic {Ai z = tt), and otherwise is something in between. 

Fig. 2 of Ref. [|| shows that the specific heat in the 
three dimensional model has a peak around T = 1.05. 
The maximum value of the plotted quantity is 24 x C ~ 
36. The result was found to be about the same for all 
three different system sizes in the figure, of which the 
smallest was 48 x 40 x 48. To compare with the specific 
heat for a single 2D layer, we have made simulations with 
L — 44 which is chosen to have close to the same area as 
48 x 40. Recall that the peak of the specific heat in the 2D 
model is a non-singular feature well separated from the 
Kosterlitz-Thouless temperature T KT ~ 0.892 jl2|. 

Our results for 24 x C for the 2D system are shown 
in Fig. ||. The similarity to the data of Hu and Tachiki 
(inset of their Fig. 2) is striking. Not only the shape 
of the curves is the same but also the position and the 
height of the maximum: at T rs 1.05 we find the value 
« 36.6 which is the same as - or possibly slightly higher 
than - the data in their Fig. 2. We therefore conclude 
that the data for rj 2 = 100 of Ref. Q is dominated by 
2D fluctuations. The absence of any sharp features in 
this data should therefore not be taken as evidence for a 
continuous transition in the 3D model. 

Why is the interplane coupling not more important? 
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An obvious but rough way to assess the importance of 
the coupling between two planes is to compare the total 
coupling with the temperature. The size of the planes 
should then be large enough if L x x L y x J z ^> T, which 
for our present values becomes 19.2 3> 1.05. This would 
lead to the conclusion that the planes are large enough to 
make the interplane coupling very significant. The reason 
why this approach fails is the presence of the frustrating 
vector potential Ai Z that varies between and 2-7T. This 
means that the interplane couplings tend to cancel each 
other out, and we believe this to be the reason for the 
very 2D-like behavior found in the model. 



line of reasoning presumes that (u 2 z ) not is strongly af- 
fected by the melting in the x direction, but there is no 
a priori reason for this assumption. Figure [3] illustrates 
that the disorder in the x direction may well give rise to 
forces that destroy the layered structure. This suggests 
problems with the stability of the smectic phase, which 
would seem to altogether invalidate the continuous melt- 
ing scenario. 



a) 
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FIG. 2. Specific heat for a two-dimensional XY model with 
L = 44. The normalization is chosen to facilitate direct com- 
parison with Fig. 2 of Ref. ||. The conclusion is that their 
data is dominated by 2D fluctuations and not should be con- 
sidered as evidence for a continuous transition in the 3D sys- 
tem. 

In the theory by Balents and Nelson [Q the transition 
takes place in two steps with an intermediate smectic 
phase. The smectic phase is characterized by order in 
one direction only, which results in only certain layers 
being occupied with vortex lines and the others being 
empty. However, to us this picture doesn't seem to be 
a convincing explanation of the experimental findings of 
Ref. H for two reasons. 

First, the experimentally found smooth onset of resis- 
tivity was only found for very small angles, 9 < 0.5°, but 
there seems to be no reason to expect that the transi- 
tion between the smectic and the liquid phases should be 
very sensitive to the precise angle of the applied magnetic 
field. The transition between the tilted smectic phase and 
the liquid is also expected to be a continuous transition, 
which only changes character to first order for large tilts 

i- 

Second, considering the thermal fluctuations of the 
vortex lattice, Balents and Nelson found that the average 
relative displacements in the x and z directions are un- 
equal with (u 2 .) /a x ^> {u 2 z ) ja? z . This was then taken to 
suggest that the Lindemann criterion for melting can be 
satisfied in the one direction only and not in the other, 
leading to a partial melting of the lattice. However, this 
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FIG. 3. The figures show a side-view of a layered super- 
conductor in the direction of the applied field. Panel a) is a 
vortex lattice and panel b) illustrates a smectic phase which is 
ordered in the z direction (only certain layers occupied) but 
disordered in the x direction. This disorder may, however, 
give rise to forces (as indicated by the arrow) that tend to 
move the vortex lines from the occupied layers and destroy 
the smectic order. 

So far we have both discussed the validity of the claim 
of a continuous transition from simulations and some 
problems with the theory of continuous melting. We now 
turn to the experimentally found smooth onset of the re- 
sistivity H suggestive of a continuous transition. We will 
argue that this observation not is inconsistent with a first 
order melting transition. 

A current through the system perpendicular to the ap- 
plied field gives rise to a force on the vortex lines perpen- 
dicular to both the current and the field. If the vortex 
lines move in response to this force the system dissipates 
energy, which is experimentally seen as a resistivity in 
the sample. In the geometry of Ref. || with the current 
in the CuO planes and the applied field also parallel to 
these planes, dissipation is related to the motion of the 
vortex lines in the z direction. 

An interesting observation is that the mechanism for 
motion of vortex lines in the z direction becomes qualita- 
tively different for perfect alignment of the field with the 
CuO layers as compared to a tilted magnetic field. This 
is illustrated in Fig. |4|. Panel a) shows the situation in 
the presence of a tilted magnetic field where the upward 
motion of a vortex line is linked to the vortices moving to 
the left in the figure. Panel b) shows the qualitatively dif- 
ferent situation for fields perfectly aligned with the CuO 
planes. In this case the motion of a vortex line has to 
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start with a 'kink' that then grows bigger until the whole 
line has crossed the plane. The crucial thing to note is 
that such a kink gives rise to a vortex anti- vortex pair in 
the x-y plane, and that this vortex pair separates as the 
vortex line gradually moves upward across the plane. 

The energy barrier to be overcome for motion of the 
vortex lines in Fig. is the pinning of individual vor- 
tices. In contrast, the energy barrier in Fig. is the kink 
energy barrier - related to the vortex anti-vortex inter- 
action in the plane - which may be quite large just above 
T c . If we assume that the transition implies a change of 
this potential barrier from infinite to quite large, the first 
order transition should not be clearly seen in the resis- 
tivity. In a 2D system this potential barrier diverges as 
T — > Tkt, and as we expect T c — > Tkt as the anisotropy 
increases, and that the behavior is 2D-like above T c , we 
believe that this assumption of a large potential barrier 
right above T c is a reasonable one at least for fairly large 
anisotropies. 




be no sharp anomaly in the specific heat. However, so 
far the experiments are inconclusive. In Ref. [ fl3| the au- 
thors found features in the specific heat both for field 
parallel and perpendicular to the planes. The work was 
later extended to several different values of the angle of 
the field, including a parallel field, and it was found that 
the sizes of these features scaled with the angle in accor- 
dance with a simple scaling relation |Q. This seems to 
point to the existence of a first order melting transition, 
but not conclusively. As the authors point out | Q , jl4| , 
the reason for the apparent first order character could 
also be that the angular resolution was not sufficient to 
access the true perfect alignment region |0| < 0.5° j|. A 
carefully designed experiment should make it possible to 
discriminate between the two different proposed pictures. 

To conclude, we have found a single first order transi- 
tion in the anisotropic 3D XY model with magnetic field 
applied parallel to the planes. The finding that very large 
system sizes are needed for obtaining the true behavior 
make us doubt the continuous transitions reported from 
simulations on smaller systems. Some weaknesses with 
the commonly accepted theory of continuous vortex lat- 
tice melting have been discussed and we suggest that the 
reason for the smooth onset of resistivity is a large poten- 
tial barrier right above T c associated with the separation 
of vortex anti-vortex pairs which is linked to the motion 
of the vortex lines in the z direction. 
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tel and Professor Petter Minnhagen. This work has 
been supported by the Swedish Natural Science Research 
Council through Contract No. E-EG 10376-312, and by 
the resources of the Swedish High Performance Comput- 
ing Center North (HPC2N). 



FIG. 4. Motion of vortices and vortex lines. These figures 
illustrate the motion of vortex lines for a) a tilted magnetic 
field and b) perfect alignment. For perfect alignment of the 
field there is a vortex anti- vortex pair in the CuO plane, and 
together with the motion of an increasingly larger portion of 
the vortex line to the upper layer goes a separation of the 
vortices in the CuO plane. 

With this mechanism for the transition the very strong 
dependence of the resistivity on the angle of the ap- 
plied field H follows naturally. Perfect alignment of 
the field corresponds to two-dimensional XY layers with 
only thermal vortex pairs and no free vortices in the low- 
temperature phase. But even a small non-zero angle of 
the field gives some additional free vortices, and the pres- 
ence of them would lead to a qualitatively different be- 
havior. 

The existence of a first order melting transition should 
be possible to verify from measurements of e.g. the spe- 
cific heat. This would be a direct experimental test of our 
proposed picture, since in the alternative scenario with a 
continuous transition into a smectic phase there should 
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